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We present a density functional study of graphene adhesion on a realistic Si02 surface taking 
into account van der Waals (vdW) interactions. The Si02 substrate is modeled at the local scale 
by using two main types of surface defects, typical for amorphous silica: the oxygen dangling bond 
and three-coordinated silicon. The results show that the nature of adhesion between graphene and 
its substrate is qualitatively dependent on the surface defect type. In particular, the interaction 
between graphene and silicon-terminated Si02 originates exclusively from the vdW interaction, 
whereas the oxygen-terminated surface provides additional ionic contribution to the binding arising 
from interfacial charge transfer (p-type doping of graphene). Strong doping contrast for the different 
surface terminations provides a mechanism for the charge inhomogeneity of graphene on amorphous 
Si02 observed in experiments. We found that independent of the considered surface morphologies, 
the typical electronic structure of graphene in the vicinity of the Dirac point remains unaltered in 
contact with the Si02 substrate, which points to the absence of the covalent interactions between 
graphene and amorphous silica. The case of hydrogen-passivated Si02 surfaces is also examined. In 
this situation, the binding with graphene is practically independent of the type of surface defects 
and arises, as expected, from the vdW interactions. Finally, the interface distances obtained are 
shown to be in good agreement with recent experimental studies. 

PACS numbers: 68.35.Np, 68.35.Dv, 73.22.Pr, 73.20.At 
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I. INTRODUCTION 

Graphene, a single layer of graphite, is subject of inten- 
sive research owing to its remarkable electronic proper- 
ties, which make this material a promising candidate for 
various electronic applicationsii"— Since graphene is one- 
atom thick, its properties are strongly influenced by the 
environment, such as, for example, atomic impurities)^ 
molecular adsorbates;^ metallic contacts^^ii and dielec- 
tric substratesi^ In respect to practical applications, the 
behavior of graphene in contact with substrates is of great 
importance. Since the discovery of graphene the role of 
substrates in its properties is actively debated, but is still 
not clearly understood. 

Silicon dioxide is the main insulating component of 
present microelectronic devices owing to its perfect di- 
electric behavior and chemical stabilityJ^ However, this 
material is not the best candidate for its use as a sub- 
strate for graphene-based devices. In comparison to the 
remarkably high carrier mobility observed in suspended 
graphene fii carrier mobility in graphene supported on 
Si02 is very limited by scattering from surface states, 
by surface phonons, and by structural deformations of 
graphene ii2r2i The deformations are attributed not to 
intrinsic corrugations of graphene; ^^'^^ which might also 
exist in Si02-supported graphene^ii but rather to the 
morphology of the Si02 surface itself, which is amorphous 
and, therefore highly disorderedJ^ It has been demon- 
strated that the rippling of graphene can be consider- 



ably suppressed by its deposition on relatively flat sur- 
faces, such as, for example, muscovite micapi^ or hexago- 
nal boron nitrideJ^ Furthermore, it has been shown that 
the adhesion of graphene on Si02 is essentially conformal, 
meaning that graphene reproduces the substrate topog- 
raphy with very high accuracy.— Apparently, interfacial 
interactions play a key role in the understanding of the 
structural behavior of graphene supported on surfaces. 
Nevertheless, the microscopic nature of the graphene ad- 
hesion on realistic Si02 has not yet been clearly estab- 
lished. 

Apart from the structural changes of graphene, a sub- 
strate can directly affect its electronic properties by 
induced charges as shown by scanning tunneling mi- 
croscopy (STM) experimentsi^ Furthermore, recent elec- 
trical field microscopy (EFM) experiments have demon- 
strated that the electrical potential of amorphous Si02 
surface is rather disordered, which suggests the existence 
of localized chargelike impurities or surface defects 3^ The 
morphology of possible surface defects, their influence on 
the properties of graphene, and the bonding mechanisms 
in the graphenc-substrate system are yet unclear from 
the theoretical point of view and require a detailed mi- 
croscopic analysis. 

So far, a number of papers reported theoretical first- 
principles investigations of graphene on Si02'^'^"— For 
modeling of the amorphous silica surface, various crys- 
talline approximations are usually used paying not much 
attention to the structural peculiarities of realistic Si02. 
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As a consequence, the results substantially differ from 
one approximation to another. For instance, the oxygen- 
terminated (001) surface of a-quartz was predicted to 
be very reactive^^^ significantly affecting the electronic 
properties of graphene, whereas the corresponding recon- 
structed surface is essentially inert and has no influence 
on the electronic structure of graphene<^ Although the 
structure of amorphous Si02 is not trivial, it has only cer- 
tain types of structural features (defects) ^Si^ and thus 
could be rationally modeled only within the frame of 
specific crystalline approximations. In addition, previ- 
ous investigations of graphene on Si02 do not explic- 
itly consider dispersive interactions, but employ only lo- 
cal density approximations (LDA) for the treatment of 
exchange-correlation effects. Although the LDA inter- 
layer distances are very close to experimental values for 
some layered materials (e.g., graphite), this approxima- 
tion is purely local and not able to describe dispersive 
interactions of nonlocal nature. As a consequence, it is 
impossible to reproduce accurately the binding energies 
and elastic properties of vdW complexes within the LDA. 
As has been shown, the dispersive interactions play a 
significant role in the adhesion of graphene on silicate 
surfacesi^ For this reason, it is desirable to use more so- 
phisticated approaches, allowing for explicit treatment of 
the vdW interactions. 

In this paper wc present a first-principles van dcr 
Waals density-functional study (vdW-DF) of graphene 
supported on amorphous Si02. Particularly, we are fo- 
cused on the interaction between graphene and its sub- 
strate as well as on the electronic structure of supported 
graphene. We model the Si02 substrate by considering 
idealized surfaces with two main surface defects, typi- 
cal for amorphous silica: the oxygen dangling bond and 
three-coordinated silicon. It is assumed that this struc- 
tural model represents a realistic Si02 surface at the local 
scale. Additionally, we also examine the effect of hydro- 
gen passivation of the surface defects considered. The 
account of the vdW interactions allowed us to estimate 
the vdW part of the graphene-Si02 interaction, which is 
found to be ranging from 23.8 to 39.6 meV/C depending 
on the particular surface type. In the case of hydrogen 
passivated surfaces as well as for the silicon-terminated 
surface, the vdW interactions play a dominant role. By 
contrast, for the oxygen-terminated surface, the vdW 
part is responsible only for half of the total binding en- 
ergy. The rest arises from the charge transfer and from 
resulting ionic interactions between graphene and the un- 
derlying surface. In addition, we show that in the ab- 
sence of interfacial covalent bonding with the substrate, 
the electronic structure of graphene remains virtually un- 
altered, which may play an important role for practical 
applications. Finally, calculated interlayer distances and 
obtained charge inhomogeneity for graphene on Si02 are 
found to be in agreement with the data of atomic force 
microscopy (AFM) and STM experiments ; respec- 
tively. This agreement justifies the structural model and 
computational approaches used in the present study. 



The paper is organized as follows. In Sec. II, we briefly 
describe the methods of the investigation. In Sec. Ill, 
we present the structural model of the Si02 surface and 
describe other structural parameters used in our study. 
In Sec. IV A, we analyze the properties of the model 
Si02 surfaces. The results on interface geometry, adhe- 
sion, and electronic properties of graphene on Si02 are 
discussed in Sees. IV B, IV C, and IV D, respectively. 
Comparison with other theoretical studies are given in 
Sec. IV E. In the last section we briefly summarize our 
results. 



II. CALCULATION DETAILS 

The results presented in this study have been ob- 
tained based on the plane-wave pseudopotential method 
as implemented in the QUANTUM-espresso simulation 
packagei^2i^ To calculate adhesion energies and prop- 
erly take into account dispersive interactions, exchange 
and correlation effects have been treated according to the 
vdW-DF method proposed by Dion et alM^ Up to now, 
this method has been extensively applied to a large vari- 
ety of compounds showing transferability across a broad 
spectrum of interactions, such as ionic, covalent, and 
vdWi^ Although the vdW-DF method docs not usually 
provide chemical accuracy, it has distinct advantages over 
the standard (semi)local approximations, such as explicit 
inclusion of the nonlocal electron correlations, leading to 
the proper description of the dispersive interactions. 

Within the vdW-DF, the exchange-correlation energy 
functional consists of several contributions, 

EUn) = Er^^^in) + E^^^{n) + Ef{n), (1) 

where the first term corresponds to the exchange part 
of the revised PBE (revPBE) functional^ is the 

correlation energy within the LDA, and i?"'(n) is the 
nonlocal correlation correction, which is calculated in the 
following way: 

Ef = dVdVn(r)^(r,r')n(r'), (2) 

where n{r) is the electronic density and (j)(r, r') is a func- 
tion incorporating many-body density response (for de- 
tails see Ref. |33 ) . It should be emphasized that we eval- 
uate the nonlocal correction [see Eq. ([2])] in a pertur- 
bative way, i.e., using only GGA-based (semilocal) elec- 
tronic density distribution. This choice looks quite rea- 
sonable since the vdW interactions are rather weak and 
thus cannot significantly change the electronic distribu- 
tion. Moreover, as has been shown on different vdW 
complexes, the effects due to lack of self-consistency are 
negligiblci^ To calculate this correction to the total en- 
ergy, we used an efficient post-processing routine noloco 
from the exciting code.— >^ 

Self-consistent treatment of the dispersion effects 
within the vdW-DF approach is considerably more de- 
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manding in terms of the computational time. As a re- 
sult, an efficient calculation of forces acting on atoms is 
also rather problematic for large systems. For this rea- 
son, we do not use straightforward geometry optimiza- 
tion of the graphene/Si02 system within the vdW-DF 
functional. Instead, we optimize only the isolated Si02 
surfaces within the gradient corrected functional in the 
revised PBE paramctrizationj^ This choice appears rea- 
sonable since the vdW interactions can be neglected for 
the energies of silica structures. Furthermore, we assume 
that the structure of graphene remains essentially unper- 
turbed upon adhesion (see Sec. IV B for discussion). In 
order to find equilibrium separations between graphene 
and its substrates, we calculated adhesion energies as a 
function of the distance between them: 

Eadh{d) = i?Gr+Surf(d) — Eq^ — Egurf , (3) 

where i?Gr+Surf(d) is the vdW-DF energy of graphene 
separated at the distance d from the surface, whereas Eq^ 
and -Esurf are the energies of nonintcracting graphene 
and the surface, respectively. The minimum of the given 
function corresponds to the adhesion energy at the equi- 
librium separation, min{£'adh((^)}=-E'adh(rfeq)- 

In our calculations, we employed an energy cutoff of 
30 Ry for the plane- wave basis and 300 Ry for the charge 
density. Self-consistent calculations of the Kohn-Sham 
equations were carried out employing the convergence 
criterion of 10~* Ry. For accurate Brillouin-zone inte- 
gration, the tetrahedron scheme^ and a (16 x 16 x 1) 
Monkhorst-Pack k-point mesh^S were used. A much finer 
mesh (48 x 48 x 1) and a Gaussian broadening of 0.01 
Ry were used for the density of states (DOS) calculations. 
We checked that further increase of computational accu- 
racy does not significantly change the results. In all cases 
under consideration, the height of the supercell was cho- 
sen to be 50 A. In order to avoid spurious interactions 
between images of the supercell in the [001] direction, we 
also used a dipolc correction!^ The convergence criterion 
for the relaxation of the Si02 surfaces was set to 0.001 
Ry/A. The positions of the lowermost layer of atoms were 
fixed. 



III. SURFACE STRUCTURES 

First-principles studies of amorphous compounds like 
Si02 are challenging due to lack of translational sym- 
metry in these systems. In many cases the surface 
structure of amorphous Si02 is modeled using various 
crystalline approximations i^^^ ^^^i^^" — Specifically, vari- 
ous faces and terminations of a-quartz and /3-cristobalite 
are frequently used to reproduce the structure of Si02. 
In this respect the structure of /3-cristobalite seems to 
be more appropriate since /3-cristobalite and amorphous 
Si02 have very similar local structures, as has been shown 
by neutron diffraction experiments However, in con- 
trast to crystalline phases of Si02, the surface structure 
of the amorphous phase is nonuniform and has different 
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FIG. 1. (Color online) Surface structures of Si02 with the 
main types of defects (side view): (a) threefold coordinated 
silicon, Q3, (b) Q3 terminated by hydrogen atom, (c) single 
oxygen dangling bond, Q\, and (d) Q\ terminated by hydro- 
gen atom. 



reactive sites, which leads to unique chemistry of this 
surface 4^ 

According to atomistic molecular dynamics 
simulations;^ there are three types of stable point 
defects on amorphous Si02 surface: (i) 03=Si», three- 
fold coordinated silicon with no dangling oxygen bonds 
(Qa), (ii) 03=Si— 0», fourfold coordinated silicon with 
one dangling oxygen bond {Q\), and (iii) threefold 
coordinated silicon with one dangling oxygen bond {Q\)- 
The notation QJ^ here implies m-coordinated silicon 
with n oxygen dangling bonds. For convenience, we 
will use this notation throughout the paper. Moreover, 
as follows from the atomistic simulations, and Q\ 
arc the most common defect sites (^--^80% of the total 
surface area), which is in agreement with other theo- 
retical investigations!^ From the experimental point 
of view these two defect types also can be considered 
as fundamental in amorphous Si02i^ An additional 
important type of chemical environment for silica is 
(5°, which corresponds to bulk silica without defects 
and can be also modeled as a passivated surface defect 
(e.g., by hydrogen). In practice, partial passivation 
of the Si02 surface can be achieved by exposure to 
molecular hydrogen (H2) even at room temperatures 
In principle, the emergence of other surface defects is 
possible, but their existence is energetically unfavorable 
and, with high probability, leads to a local reconstruction 
of the surface. Interestingly, the reconstructed surface 
of crystalline a-quartz phase is predicted to be fully 
saturated owing to the formation of six-mcmbcred rings 
at the surfacCi ^^'^^ This is. however, a consequence of 
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FIG. 2. (Color online) Schematic representation of graphene 
supported on a Si02 surface with different defects on the top 
(according to the model employing in this work). 



the long-range order of crystalline surfaces and is not 
likely to occur for amorphous surfaces. 

Unsaturated surfaces arc usually reactive and eas- 
ily adsorb impurities from the environment to passi- 
vate themselves. Such a passivation can be, in princi- 
ple, avoided under ultrahigh vacuum (UHV) conditions. 
However, as has been previously demonstrated by EFM 
experiments, the amorphous Si02 surface shows rela- 
tively large fluctuations of electrical potential at ambi- 
ent conditions)^ This fluctuations suggest the existence 
of uncompensated surface charges even without UHV, 
which can be associated with unsaturated surface defects. 

We believe that the chemical reactivity of realistic 
amorphous silica can be defined by a certain combination 
between the reactivities of local surface regions. This re- 
gions can be modeled by idealized surfaces, each of which 
contains only one defect type among those discussed 
above. In turn, the properties of graphene at a large scale 
(reahstic Si02) are expected to be composed of its prop- 
erties at the local scale (idealized Si02). Besides compu- 
tational advantages within the first-principles framework, 
this approach allows us to investigate the role of different 
surface defects separately. In our work, for modeling the 
Si02 surface, we consider two different /3-cristobalite sur- 
faces with morphologies corresponding to the Q3 and Q\ 
defect sites. We show relaxed structures of these surfaces 
in Figs. [TJa) and[TJc). As we already mentioned, these 
two defects are the most numerous types in amorphous 
Si02 according to the previous studies.— 1^ We do not 
consider the Q3 defect mainly because it would have been 
necessary to consider a (100) surface of /3-cristobalite in- 
stead of a (111) surface that is suitable for the modeling 
of and Q^-defective surfaces. However, the unit cell of 
a (100) /3-cristobalite surface is not commensurate with 
the unit cell of graphene in lateral directions. In this 
case, the lattice constant mismatch is more than 10% 
and the corresponding adjustment of the Si02 unit cell 
would considerably affect the properties of the surface. 



Since the Ql defect contains one oxygen dangling bond 
and one three-coordinated silicon, it is reasonable to sup- 
pose that this defect type may be approximated by two 
separated defects, Ql and Q^. In addition to unsaturated 
surfaces, we investigate the effect of hydrogen passivation 
on adhesive properties of Si02 surfaces [see Figs.fljb) and 
[IJd)]. In Fig. [31 we show schematically the structure of 
graphene supported on Si02 with various defects on its 
surface. 

In all cases, the supercell used in our study contains 
three layers of dcfectless (Si04)"'^ tetrahedra and one 
topmost defective layer. Oxygen atoms in the lower- 
most layer are saturated by hydrogen atoms. All atoms 
except the lowermost layer of hydrogen were relaxed. 
We used lateral (2 x 2)a graphene unit-cell parameters 
for the Si02-graphene supercell, which corresponds to a 
^3% laterally compressed (111) surface unit cell of /3- 
cristobalitei^ We believe this choice is justified since sil- 
icates are known to have a relative low bulk moduli, i.e., 
they can be compressed quite easily<^ 

As for graphene, it has a two-dimensional honeycomb 
lattice of sp^-bonded carbon atoms. Although the real 
structure of graphene is corrugated, the characteristic 
length of corresponding ripples is around 100 A^^^ 
which is much larger than the typical length of super- 
cclls used in first-principles calculations. For this reason, 
we do not consider this phenomenon in our work. We 
used a lattice constant of graphene equal to a = 2.459 A 
in accordance with the experimentally obtained value for 
graphite at low temperaturesi^ 



IV. RESULTS AND DISCUSSION 

A. Bare Si02 surfaces 

In Fig. 131 we show densities of electronic states for the 
Si02 surfaces under investigation. One can see that the 
H-passivated surfaces [see Figs. M)^) and ^d)] are fully 
saturated and there are no resonances within the band 
gap, indicating the absence of localized states at the sur- 
face. In these cases, the DOS exhibits a band gap of 
~5.7 eV, which is comparable with the band gap of bulk 
/3-cristobalite calculated using the PEE functional (5.84 
eV)i^ These values, however, underestimate the experi- 
mental gap of amorphous Si02 films (9.3 eV)r^ which is 
related to the well-known problem of the standard den- 
sity functional theory (DFT) functional (LDA/GGA), 
such as inability to accurately describe electronic spec- 
tra for compounds with strong electron localization. The 
band gap of silica can be considerably improved using hy- 
brid functionals (e.g., PBEO)^ or GW approximation!^ 
In our work, however, we are not focused on a precise 
quantitative description of the electronic properties of 
Si02. We believe that the use of more appropriate ap- 
proaches will not qualitatively affect the results obtained 
in this study. 

In contrast to H-passivatcd Si02, defective surfaces are 
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FIG. 3. (Color online) Densities of states for Si02 surfaces 
projected onto Si (blue line) and O (green line) atoms; (a) 
and (c) surfaces covered by the Q3 and Qi defects, respec- 
tively; (b) and (d) corresponding hydrogen-passivated sur- 
faces. Oxygen DOS is reduced by a factor of five for clarity. 
The vertical line corresponds to the Fermi level. Zero energy 
coincides with the lower energy of the conduction band. 



TABLE I. Work functions (WF) and electron affinities (EA) 
calculated for different Si02 surfaces (see Sec. Ill for details) 
and graphene. 





Ql 


Ql + H 


Ql 


Ql + H 


graphene 


WF, eV 


5.6 


8.0 


9.3 


7.2 


4.2 


EA, eV 


4.1 


2.3 


8.6 


1.6 


4.2 



unsaturated and, therefore the DOS exhibits resonances 
within the band gap arising from the localized surface 
states. For each surface, there are two peaks correspond- 
ing to occupied and unoccupied states of opposite spin. 
In the case of surface with undercoordinated silicon (the 
Qg defect), the midgap states originate from the Si(3s,3p) 
and 0{2p) orbitals, whereas for the Q\ defect the reso- 
nance is slightly more intense and localized mainly at the 
0{2p) orbital. The presence of localized surface states 
allows one to expect chemisorption of graphene on a sur- 
face. Indeed, with respect to graphene supported on 
Si02, the defect sites may be effectively considered as 
monovalent impurities. In turn, such impurities may bind 
ionically or covalently with graphenei^ 

In order to estimate chemical activity of the surfaces 
in respect to graphene, we calculate the work func- 
tions (WF) and electron affinities (EA) as (Evac-Ep) and 
(^'vac-^'cond), respectively, where E^^c is the vacuum level 
of the electrostatic potential, Ep is the Fermi level, and 
-Bcond is the lowest energy of the conduction band. We 
stress, however, that the accuracy of the WF and EA 
is limited within the employed computational method, 
similar to the band-gap problem (see above). Neverthe- 
less, these quantities arc useful to establish general trends 



in the chemisorption of graphene on Si02. In Table HI 
we summarize calculated WF and EA for Si02 surfaces 
and graphene. For graphene, WF is equal to EA due to 
the absence of a band gap. Comparing WF and EA for 
graphene and H-passivated surfaces, one can sec that the 
surface WF is larger than the EA of graphene (and vice 
versa). The same situation occurs for graphene and the 
surface with the Q3 defects. This implies that the charge 
transfer between these surfaces is not energetically favor- 
able. On the other hand, the WF of graphene is lower 
than the EA of the surface with oxygen dangling bond 
(Ql), which suggests electronic transfer from graphene 
toward this surface (p-type doping of graphene) . 



B. Geometry of the graphene-Si02 interface 

As we already mentioned in Sec. II, we calculate equi- 
librium graphene-Si02 distances assuming that the struc- 
ture of graphene is unperturbed upon contact with the 
substrate. However, this does not reflect the structure of 
graphene on Si02 on a large scale. As follows from ex- 
perimental studies, the structure of graphene supported 
on amorphous Si02 is highly corrugated owing to high- 
fidelity conformation! ^^-^^1^° The assumed invariability 
of graphene is resulting from the employed theoretical 
approach. Indeed, the supercell used in our study is 
too small in lateral directions to reveal corrugations of 
graphene caused by the irregularity of the substrate. 
Moreover, imposed boundary conditions along with the 
strong sp^ bonding of carbon atoms artificially restrict 
the flexibility of graphene. Previous results show that 
for such small supercells the deformation of graphene is 
really negligible!^ Thus we note that our results reflect 
only the local behavior of graphene on Si02 . Therefore 
the interlayer distance is the main structural parameter 
in our model. In Table |lll we summarize these parame- 
ters for different Si02 surfaces. 

Equilibrium distances between graphene and Si02 are 
considerably different for the surfaces covered by the 
Q3 and Q\ defects. This can be easily accounted for 
by noticing that the difference in the distance correlates 
with the difference in the adhesion energy. Namely, the 
smaller distance corresponds to the larger i^adh and vice 
versa. The situation for hydrogen-passivated defects is 
different. In this case, although the binding energy is 
smaller relative to the binding for the defect, the in- 
terlayer distance is also smaller. This is caused by the 
fact that the orbitals associated with hydrogen have a 
small spread, i.e., they are not much spatially extended 
compared to the oxygen or silicon orbitals. As a conse- 
quence, the Pauli repulsion between graphene and the hy- 
drogen layer decreases more rapidly with distance, which 
allows graphene to come closer to hydrogen. 

The interface distances obtained between graphene 
and Si02 can be compared with experimental height his- 
tograms acquired across the graphene-Si02 boundary us- 
ing the AFM technique According to the experimental 
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TABLE II. Calculated interlayer distances, adhesion energies, 
Fermi level shifts, and WFs (EAs) for graphene supported on 
different Si02 surfaces. Interlayer distances are calculated 
as a difference between the z coordinate of carbon atoms in 
graphene and the topmost atomic layer of the substrate. The 
Fermi level shift implies the difference between the Dirac point 
and the Fermi level of supported graphene. The values are 
given for the most stable graphene geometry. 





Q% 


Ql + H 


Qi 


Q\ + H 




3.6 


2.7 


2.7 


2.6 


£adh, meV/C 


-31.6 


-26.5 


-11.1 


-24.8 


vdW part of -Eadi^l 


98% 


96% 


51% 


96% 


I^Ef, eV 


0.0 


0.0 


-1.0 


0.0 


WF (EA), eV 


4.2 


4.2 


5.2 


4.2 



Contribution of vdW interaction calculated as a difEerence 
between adhesion energies in the presence of nonlocal 
correlations and without it. 



data, the average thickness of a single graphene layer de- 
posited on Si02 under UHV conditions can be estimated 
as 4.2 A. This quantity, however, does not reflect the 
actual distance of the interface since it directly depends 
on the spatial extension of the orbitals of both surfaces. 
In fact, the tt orbitals of graphene are known to have a 
rather large spread (see, e.g., Ref. [60[) and, therefore, the 
thickness of graphene supported on a surface with more 
localized orbitals should be larger than the distance in the 
theoretical sense, i.e., as a difference between atomic co- 
ordinates. On the other hand, if the thickness of the top- 
most atomic layer is much smaller than the thickness of 
supported graphene, the interlayer distance corresponds 
to half of the graphene thickness. It is difficult to estab- 
lish a more precise relationship between these quantities 
because of a number of ambiguities in the determination 
of the orbital spread and its quantitative influence on ex- 
perimental results. According to this consideration, the 
average experimental distance should satisfy the follow- 
ing uncertainty relation: hjl < dc-^p < h, where h is the 
average experimental thickness. Thus, for the mentioned 
experiment^ the distance dexp is bounded within 2.1- 
-4.2 A. One can see that obtained theoretical distances 
calculated according to the given structural model of the 
Si02 surface (see Table HI]) fall entirely within this range. 
Taking the uncertainty into account, we come to a good 
agreement between results obtained and experimental es- 
timations. It should be noted that the obtained distances 
are rather sensitive to the presence of the vdW interac- 
tions, especially if they dominate. In fact, the distances 
calculated without the nonlocal term in the exchange- 
correlation functional [see Eq. are significantly higher 
and therefore do not lead to an agreement with the ex- 
periment. 



C. Graphene adhesion on Si02 



As follows from Table HIl the main contribution (96%) 
to adhesion of graphene on H-passivated Si02 surfaces 
comes from the vdW interaction. Adhesion energies for 
Q'^ + H and Q\-\- H cases are very close to each other, 
which indicates that the vdW part of the adhesion en- 
ergy is only slightly dependent on the particular surface 
type. The situation for unsaturated surfaces is partially 
different. In the case of a surface with undercoordinated 
silicon (Qs), the vdW contribution still dominates (98%) 
leading to the adhesion energy that is slightly larger than 
for the H-passivated case, whereas the presence of dan- 
gling oxygen bonds on the surface (Q4) results in a much 
stronger adhesion of only 51% of the vdW interaction. 
The latter implies that graphene definitely chemisorbs 
on the surface with Q\ defect, which is the most reac- 
tive among the others. The values in Table [ll] are given 
for the most stable configuration of graphene relative to 
the considered surfaces. Namely, for all considered sur- 
face defects, this configuration corresponds to the situ- 
ation where point defects are situated below a hollow 
(hexagon center) graphene site. If the defect is located 
below a bridge (center of the C-C bond) or top site (C 
atom), the maximum difference in adsorption energy is 
reached for the oxygen terminated surface (^^5 meV/C). 
For three other surfaces considered, this difference is not 
exceeding 1 meV/C. It is worth mentioning that irrespec- 
tive of the surface type, values obtained for the adhe- 
sion energies are sufficiently exceeding the elastic energy 
stored in graphene (^^0.8-2.6 meV/C)] ^^i^° This means 
that the interaction energy between graphene and Si02 is 
large enough to overcome the energy needed for graphene 
to follow the Si02 morphology. Therefore, our results 
confirm previous estimations regarding the possibility of 
graphene-Si02 conformation. 

Very recently, the first direct measurements of the 
adhesion energy of 1-5 layer graphene to Si02 were 
reportedi^ For monolayer graphene the value of 0.45 ± 
0.02 J/m^ (~ 74 ± 3 meV/C) was obtained. As can be 
seen, the experimental adhesion energy is very close to 
our theoretical estimation (77.7 meV/C) for the oxygen- 
terminated surface {Q\ defect), but significantly higher 
than for the other surface types. This might imply that 
the Q\ defect is the most common defect at the Si02 
surface in the experiment. However, in view of the fact 
that the Q\ defect can be considered as a charge impu- 
rity (as we show below), this assumption seems to be 
inconsistent with the reported expectations that charge 
impurities should not significantly affect the adhesion be- 
cause of their relatively low density in Si02. We believe 
that in order to resolve this contradiction, a more reli- 
able large-scale model for the Si02 structure is required, 
where concentration and distribution of the (sub) surface 
defects are explicitly taken into account. We leave this 
problem open for further studies. 



D. Electronic structure and charge transfer in 
supported graphene 

In Fig.m we show DOS for graphene on Si02. As there 
is no chemisorption between graphene and H-passivated 
Si02 surfaces, the electronic structure of these systems 
is not changing upon adhesion [see Figs. SJ^b) andSlJd)]. 
In this situation, the Dirac point of graphene coincides 
with the Fermi level and lies in the center of the sub- 
strate band gap. The same is true for the Q^-defective 
surface with the exception that the DOS of graphene is 
slightly distorted in the vicinity of impurity resonances. 
This may indicate a slight orbital hybridization between 
graphene and the surface. However, these distortions are 
small enough to produce any noticeable changes in the 
adhesion energy. As we have shown before, in the case 
of graphene adherent to the Q\ surface, there should be 
additional sources of the interaction apart from the vdW 
one. Namely, as follows from Table U p-type doping of 
graphene is expected since the EA of the Q\ surface is 
larger than the WF of isolated graphene. This is also 
clear from the DOS [see Fig.lUc)] that exhibits the Fermi 
level shift toward the low-energy range relative to the 
Dirac point. Such a shift means that graphene works 
as a donor on the Q4-defective Si02 surface. As can be 
also seen in Fig.Ul^c), the Fermi level lies at the impurity 
resonance suggesting that the electrons from graphene 
are transferred to the p-orbital of oxygen, which becomes 
partially occupied. Lowdin charge analysis^ shows the 
charge doping of 0.29 e~/cell. 

The Fermi level shift implies the existence of induced 
charges in local regions of supported graphene, which 
can be considered as Coulomb impurities. This finding 
is in agreement with experimental STM studies, where 
electron-density inhomogeneities have been clearly ob- 
served in graphene flakes placed on a Si02 substratci^ 
Raman spectroscopy studies also provide some evidences 
for the p-typc doping of graphene on Si02'^ Our results 
allow one to identify more specifically the origin of the 
charge inhomogeneities and associate them with the Qi 
surface defect. The presence of these impurities plays a 
significant role for charge-carrier scattering and may af- 
fect the electron mobility of a sample. This is especially 
important for graphene supported on a Si02 substrate, 
where mobility is considerably limited in comparison to 
suspended graphenci^ Although there are a number of 
sources restricting the mobility of graphene, charged im- 
purities seem to be one of the dominating factors, at least 
at low temperatures.— It is interesting to note that the 
overall mobility of graphene depends on the distribution 
of the charged impurities on its surfacei^i^ In partic- 
ular, the contribution of impurities to the resistivity is 
maximal for their homogeneous distribution and can be 
strongly suppressed by clusterization. In the case under 
investigation, these effects depend on the particular Si02 
morphology. 

Apart from the Fermi level shift indicating the charge 
transfer from graphene to the Q\ surface, there is a dis- 
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FIG. 4. (Color online) Densities of states of Si02 surfaces in 
contact with graphene sheet projected onto C (black line). Si 
(blue line), and O (green line) atoms. The Si02 substrates are 
given in the same order as in Fig. [S] Oxygen DOS is reduced 
by a factor of five. The vertical line corresponds to the Fermi 
level. Zero energy coincides with the Dirac point of graphene. 



tortion of the graphene DOS within the energy range 
of —3 to 1 eV relative to the Dirac point. It is worth 
mentioning that no distortions occur in the vicinity of 
the Dirac point of graphene. The alteration of the DOS 
reflects some hybridization between the tt orbitals of 
graphene and p orbitals of oxygen from the defective layer 
of Si02. Since the oxygen atom has only one unsaturated 
electron in the Ql defect, this defect site can be consid- 
ered as a monovalent impurity in respect to graphene. 
According to Ref. 0, monovalent adsorbates on graphene 
can be divided into two separate groups regarding the 
bonding mechanism: ionically and covalently bonded im- 
purities. The main qualitative difference between these 
two groups is that the second type of impurities induces a 
band gap and corresponding midgap states at the Fermi 
level, whereas the first type gives rise only to impurity 
states lying below (or above) the Dirac point. Impor- 
tantly, the typical band structure of graphene is pre- 
served for ionic group of impurities. In accordance with 
this classification, the Q\ defect site corresponds to an 
ionic impurity since it does not disturb significantly the 
band structure of graphene. We note that no qualita- 
tive changes in terms of the total DOS have been found 
among the different surface geometries of graphene. 

Besides the covalcnt bonding, there is another mech- 
anism of the band-gap opening in graphene, resulting 
from the symmetry breaking. If a binding partner inter- 
acts with one sublattice of graphene more strongly than 
with another sublattice, the graphene sublattice symme- 
try is broken, resulting in the band-gap opening. One of 
the prominent examples exhibiting this effect is graphene 
supported on the surface of hexagonal boron nitride.— 
For graphene supported on the Si02 surfaces, we also 
found this effect, however, in all cases considered the in- 
teraction between graphene and its substrate is relatively 
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weak and the band gap is not exceeding 15 meV. Since 
this value is less than the thermal energy at room tem- 
perature (~25 mcV), we believe that this gap can be 
considered as negligible. 

In Table |lTl we also give WFs of supported graphene. 
Since the substrate does not induce a band gap in 
graphene, its WFs coincide with EAs, as for isolated 
graphene. One can see that for all considered substrates, 
with the exception of the oxygen-terminated Si02 sur- 
face, the WF of supported graphene is equal to the WF 
of isolated graphene. This is a consequence of the ab- 
sence of the chemical bonding between graphene and 
these surfaces. The WF (EA) of graphene on the sur- 
face with Ql defects is larger by the value of the Fermi 
level shift caused by the charge transfer. On a large scale, 
this means that the graphene sheet deposited on amor- 
phous Si02 should possess unequal reactivity for differ- 
ent surface regions. This property might find application 
in practice as a novel graphene-based method for micro- 
scopic surface studies, along with existing onesi^ 



E. Comparison with other theoretical studies 

Our results, in particular concerning the oxygen- 
terminated surface {Ql defect), might appear qualita- 
tively inconsistent with previous DFT investigations of 
graphene on Si02 surfaces i^i^i^^ In these works, an 
evident covalent bonding has been observed between 
graphene and Si02 along with significant structural re- 
construction both of graphene and the surface. Moreover, 
according to these investigations the electronic structure 
of graphene undergoes significant modification upon con- 
tact with Si02. Considerable band-gap opening is also 
reported. It should be noted that the results mentioned 
have been obtained using the structure of a-quartz for 
modeling the Si02 surface. This means that a cut of this 
structure in the [001] direction corresponds to a surface 
with two oxygen dangling bonds per unit cell or, in other 
words, to a (54-defective surface. In contrast, our results 
for the oxygen-terminated surface have been obtained us- 
ing the surface with only one oxygen dangling bond {Q\ 
defect) per unit cell. Therefore, the mentioned inconsis- 
tency is not surprising since the Ql defect is apparently 
more reactive than the Q\. This is also consistent with 
Ref. where no strong binding with graphene has been 
found in the case of Q^-like surface termination. As we 
already stated in Sec. Ill, the existence of free-standing 
surface configuration like Q^ is not energetically favorable 
and most probably would lead to a local reconstruction. 
Interestingly, if graphene is supported on a fully recon- 
structed Si02 surface,— neither covalent nor ionic inter- 
facial interactions occur, which means that graphene ph- 
ysisorbs on the surface. The qualitative difference in the 
properties of graphene on various Si02 surfaces accentu- 
ates the importance of the realistic surface morphology 
in the determination of graphene properties. The results 
of electronic structure calculations for the hydrogenated 



surfaces considered in our work are in qualitative agree- 
ment with the mentioned investigations. Additionally, 
we note that the previous results were mainly obtained 
using the LDA functionals and, therefore, some quanti- 
tative discrepancies (e.g., in binding energies or interface 
distances) are legitimate. 



V. CONCLUSION 

We have investigated the nature of interfacial interac- 
tions between graphene and amorphous Si02 from first 
principles. We modeled Si02 surface assuming the ex- 
istence of certain types of surface defects that are rel- 
evant for amorphous silica surface. Particularly, ideal- 
ized surfaces with the oxygen dangling bond and three- 
coordinated silicon defects are considered, as well as the 
case of their hydrogen passivation. We found that the 
vdW interaction plays a dominant role in the adhesion of 
graphene on silicon- and hydrogen-terminated surfaces. 
This confirms the importance of the dispersive interac- 
tions for graphene supported on surfaces and also justifies 
the use of the nonlocal vdW correction. For the surface 
with dangling oxygen bonds, we found, apart from the 
vdW, an additional ionic contribution to the graphene- 
Si02 binding resulting from the electron transfer from 
graphene to the more electronegative Si02 surface. This 
contribution is responsible for half of the total binding 
energy. 

We found that the DOS of supported graphene remains 
almost unperturbed retaining its typical Dirac shape. 
As a consequence, there is no covalent bonding between 
graphene and the surfaces considered. The only change 
in the electronic structure of graphene concerns the Fermi 
level shift for the oxygen-terminated substrate. This 
shift, caused by the Q\ defect, implies the existence of 
charged impurities and may be partially responsible for 
the known transport anomalies of graphene on Si02. 

Our results show that different surface terminations 
lead to the variations of the distances between graphene 
and Si02. Unfortunately, these variations cannot be di- 
rectly related to experimental observations due to the 
complexity of the realistic Si02 morphology. Neverthe- 
less, obtained interface distances are comparable with 
the average thickness of supported graphene observed 
in AFM experiments. This thickness allows us to es- 
timate the average graphene-Si02 distance. Taking into 
account ambiguities in the accurate determination of this 
distance, our theoretical results agree well with experi- 
mental data. This agreement shows the rationality of the 
structural model for amorphous Si02 used in our study. 
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